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Abstract. In this proceedings, we report measurements of the third harmonic coefficient of 
the azimuthal anisotropy, vs, known as triangular flow. The analysis is for charged particles near 
midrapidity in Au+Au collisions at ^sjvjv = 200 GeV, based on data from the STAR experiment 
at the Relativistic Heavy Ion Collider. Triangular flow as a function of centrality, pseudorapidity 
and transverse momentum are reported using various methods, including a study of the signal 
for particle pairs as a function of their pseudorapidity separation. Results are compared with 
other experiments and model predictions. 



1. Introduction 

The study of azimuthal anisotropy, based on Fourier coefficients, is widely recognized as an 
important tool to probe the hot, dense matter created in heavy ion collisions [1]. 

= (cosn((/)- ^-i?)) (1) 

where <j) denotes the azimuthal angle of the outgoing particles, is the orientation of the 
reaction plane and n denotes the harmonic. 

The first harmonic coefficient vi, called directed flow, and the second harmonic coefficient V2, 
called elliptic flow, have been extensively studied both experimentally [2] and theoretically, while 
higher even-order harmonics have also garnered some attention [3J. In contrast, odd harmonics 
of order three and above were overlooked until recently, because in a picture with smooth initial 
overlap geometry, it had been assumed that higher-order odd harmonics are required to be zero 
by symmetry. Now it is realized that event-by-event fluctuations break this symmetry [H El |6] . 
As a consequence, higher-order odd harmonics carry valuable information about the initial state 
of the colliding system [3 El [9l [lOl [TTl [121 IS] • The third harmonic coefficient — related to a 
collective motion known as triangular flow — is thus a new tool to study initial-state fluctuations, 
and the subsequent evolution of the collision system. Theoretical studies suggest that V3 is more 
sensitive to viscous effects than V2 because finer details of the higher harmonics are smoothed 
more by viscosity. We present measurements of ^3 vs. the pseudorapidity separation Ar] = 
rji — rjj between the two particles, fit with narrow and wide Gaussians, and present results for 
the wide Gaussian as a function of pT and centrality. 



2. Methods and Analysis 

Azimuthal correlations not related to the initial geometry and reaction plane orientation 
which can arise from resonances, jets, strings, quantum statistics effects, final state 




Figure 1. v'^{2,Ar]} vs. Arj for charged hadrons within two centrality intervals in 200 GeV 
Au+Au colhsions. Data are fit with narrow and wide Gaussians. Like sign (LS), unlike sign 
(US), and charge independent (CI) cases are shown. 



interactions (particularly Coulomb effects), and momentum conservation are refereed as non- 
flow correlations. Two-particle correlations which are sensitive to non-flow correlations, and 
multiparticle correlations which are less sensitive to non-flow correlations, along with event 
plane methods have been extensively used in elliptic flow measurements for the past couple of 
decades. These methods can be extended to measure triangular flow. The event plane method 
and the two-particle cumulant method are used in the present study. 

2.1. Event Plane Methods 

In the standard event plane method [1] for ^3, we reconstruct a third-harmonic event plane from 
Time Projection Chamber (TPC {\r]\ < 1.0)) tracks and also from Forward TPC( 2.5 < \rj\ < 4.0) 
tracks. The event plane vector Qn and the event plane angle ^'3 from the third harmonic of the 
particle azimuthal distribution are defined by the equations 

Q3Cos(3^'3) = = '^WiCos{34>i), (2) 
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Q3sin(3^'3) = Qsy = ^ Wi sm{3(pi) , (3) 
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Figure 2. Width and amplitude of the wide Gaussian as a function of centrahty for charge- 
independent (CI) and like-sign (LS) particles. The plotted errors are statistical. 



tan-i(^)/3, 



(4) 



In event plane calculations, tracks have a weighting factor w = pT 'vd. units of GeV/c for px < 2 
GeV /c, and w = 2 for pt ^ GeV/c. This helps to maximize the event plane resolution and 
also helps to reduce the contribution from jets in event plane calculations. Although the STAR 
detector has good azimuthal symmetry, small acceptance effects in the calculation of the event 
plane azimuth were removed by the method of shifting [T3|. When using the TPC event plane, 
we used the ry sub-event method with an additional rj gap of ± 0.05. 
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This avoids self-correlations because the particles and the event plane are in opposite 
hemispheres. When using the FTPCs, we used the full event plane from both FTPCs [1]. 
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This introduces a large rj gap between the particles and the event plane. Since there is no overlap 
between the coverage of the TPC and FTPCs, there is no possibility of self-correlation when 
using the FTPC event plane. 



2.2. Two-Particle Correlation Method 

We study (cos3((/>j — vs Ar/ for two particles with indices i and j to understand the Ar/ 

dependence of the triangular flow signal and to distinguish among different sources of non-flow 
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Figure 3. The upper panel shows the width, and the lower panel shows the amplitude of the 
wide Gaussian, as a function of transverse momentum, for the most central (0-5%) and for 
mid-central (30-40%) collisions. The plotted errors are statistical. 



correlations. This distribution of (cos3((/>j — (f)i))i^j vs. A?] can be well described by wide and 
narrow Gaussian peaks as shown in Fig. [1] for two centrality intervals. The narrow Gaussian 
is identified with short range non-flow correlations like Bose-Einstein correlations, resonance 
decay. Coulomb interactions, and effects from track merging. The narrow peak disappears above 
Pt > 0.8 GeV/c so it is unlikely to be from jet correlations. The wide Gaussian represented by 
v^{2, A?]} is the signal of interest and its fit parameters are used to calculate i'3{2} as a function 
of centrality and transverse momentum. For ^3 integrated over pT and rj we have 

^ JM{2 A„).y.(A,) 

where W equals dN/d{Ar]) when weighted with the number of particle pairs, or is set to 1 for 
unit weight. The quantity V3{2}{pt) can be obtained from the scalar product jl5j relation 

v,{2}{pt) = (^°^3('/'^--'^-))^^.-, (8) 



v/(^i{2}) 



where the j particle is selected from the pT bin of interest. 
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Figure 4. The third harmonic coefficient as a function of centrality from different methods of 
measurement for Au+Au collisions at s/sJiN = 200 GeV, with track selections 0.15 < pT < 2.0 
GeV/c and -1.0 < r] < 1.0. 



3. Results and Discussion 

Figure [2] shows the width and amplitude of the wide Gaussian as a function of centrality. We 
observe that the width peaks in mid-central collisions whereas the amplitude increases from 
central to peripheral collisions. Figure E] shows the width and amplitude of the wide Gaussian 
as a function of pT for centralities 0-5% and 30-40%. Above pT > 0.8 GeV/c, the distribution 
can be described by a single wide Gaussian. The amplitude increases with pT and then saturates 
around = 3 GeV/c. The px dependence of the width seems to depend on centrality, with 
the 0-5% most central data showing first an increase and then a gradual decrease, while at 30- 
40% centrality, the data appear to gradually decrease as a function of pT- Figure S] shows the 
centrality dependence for p^-integrated from several different analyses: f3{2} from Eq. [7]and 
Fig. [T]for the wide Gaussian, t;3{TPC} measured with third harmonic event plane reconstructed 
in the TPC, two-particle cumulants with a minimum pseudorapidity separation between particles 
of one unit, and fslFTPC} where is measured relative to the third harmonic event plane 
reconstructed in the FTPCs. 

Clearly the various analysis methods for differ greatly in Fig. [H For the wide Gaussian and 
the TPC event plane, the results are similar, suggesting that the narrow Gaussian is not very 
important. When a large Arj is required, the results decrease, especially for peripheral collisions. 
The variation between results shown in Fig. S] is caused by the Arj dependence as shown in Fig.[5j 
Non-flow effects, which are assumed to be largest when Arj is small, are expected to contribute 
less to the methods sampling a larger (|r/|). This hypothesis is consistent with previous studies 
of elliptic flow, based on two-particle correlations, where the corresponding wide Gaussian was 
ascribed to mini-jet correlations [27]. The decrease with Arj of i'^{2} may arise from decreasing 
non-flow or from decreasing initial-state density fluctuations; it is argued that those fluctuations 
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Figure 5. The third harmonic coefficient as a function of mean A?]. The two points at 
(Ar/) = 0.63 are from the methods using the TPC with \r]\ < 1. The points at 1.33 are from the 
cumulant results with |A?7| > 1. The points at 3.21 are from correlations with the FTPC event 
plane. 

decrease with increasing Arj separation |25j. Since depends so strongly on the r] separation of 
the particles, one must always quote Arj for each ^3 measurement, and compare results to models 
with the same Arj. If V3{2} is related to the initial eccentricity fluctuations, then the reduction 
of v^{2} at large Arj would presumably require a decrease of the initial state fluctuations at 
large rapidity separations. Recent work has found such a decoherence effect with a hadron and 
parton cascade model [25] . 

The non-flow contributions due to short-range correlations are effectively suppressed by using 
the wide Gaussian or by the rj-gap method. We studied the influence of different pseudorapidity 
gaps. However, Fig. [5]shows that we did not find a pattern where V3 stabilized at a constant value 
for large Aij. We compared like and unlike charge-sign combinations, since they have different 
contributions from resonance decays, fluctuations, and final state interactions. The large 
discrepancy between the methods may have its origin in the Ar] dependence of the fluctuations, 
which seem to decrease with increasing Arj separation. Thus it is not clear if one should 
extrapolate to large Ar] to avoid non-flow, or to small Ar] to measure all the fluctuations [25j. 
More theoretical input is necessary to completely understand the Ar] dependence of this signal. 
Figure [6] shows the r] dependence of V3 using event plane methods. For particles in the TPC 
using the r] sub-event method, is somewhat peaked at mid-rapidity. With the event plane in 
the FTPCs, there is a large r] gap and is flat at all centralities. This implies that acceptance 
effects at the edges of the TPC are not important. Thus, even though a large Ar] means that 
one of the particles is probably at large r/, this apparently does not strongly influence the Ar] 
dependence. 

In Fig. [71 STAR results with the event plane in the TPC are very similar to those of 
PHENIX [18J. This is surprising because the mean r] of their RXN detector is larger than 
in the case of the sub-events in the STAR TPC. However, the STAR FTPC results are lower 
than the PHENIX results. This is expected, because the mean Ar] is considerably larger in the 
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Figure 6. The third harmonic coefficient as a function of rj for different centralities for Au+Au 
colUsions at ^snn = 200 GeV, with track selection in the TPC of 0.15 < px < 2.0 GeV/c 
and \r]\ < 1.0. Results are shown for the event plane constructed both in the TPC and in the 
FTPCs. The horizontal lines are fits to the FTPC results. 



FTPC than in the PHENIX RXN detector. 

Third harmonic flow coefficient has been studied in event-by-event ideal hydrodynamics [T^] 
with MC Glauber initial conditions. These authors concluded that instead of averaged initial 
conditions, event-by-event calculations are necessary to compare with experimental data. The 
first prediction of V3 with viscous hydro was Ref. [TO]. In Fig. [51 V2{TPC} and V3{TPC} as 
a function of transverse momentum obtained with TPC sub event plane method are compared 
with several models for 0-5% and 30-40% central collisions. The experimental results for the 
TPC sub-event plane method which may have small non flow contribution are shown because 
they eliminate the short-range correlations but yet have small |Ary| like the theory calculation. 
The specific models are the viscous hydrodynamic model of Ref. |12] . where the ratio of viscosity 
to entropy is r]/s = 0.08 and 0.16, and the AMPT model [13]. AMPT model results are for 
320 < Npart < 360 and 80 < Npart < 120 from Ref [4] compared with 0-5% and 30-40% central 
data . Predictions for V3 from Parton-Hadron String Dynamics [20j at 30-40% centrality for 
\r]\ < 0.5 have been made by the sub-event method with the event planes at 1.0 < \rj\ < 4.0, and 
are also plotted in Fig. [8] (lower right panel) . 

Elliptic Q.ow{v2{T PC}) results are mostly described by ideal hydrodynamics in the case of 
the most central collisions, and by r//s = 0.08 in the case of mid-central collisions. We find that 
the third harmonic coefficient( V3{TPC}) results are also described by this model with similar 
viscosities. The PHSD model also agrees with data. 
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Figure 7. The third harmonic coefficient as a function of pT for both event plane methods 
for different centrahties compared to the results from the PHENIX experiment [18j . PHENIX 
requires \r]\ < 0.35 while STAR requires \r]\ < 1.0. In the case of STAR results from the TPC, 
the mean Aij was 0.63, while in the case of the FTPC event plane, the average Aij was 3.21. 
The PHENIX results used the event plane from their RXN detector at an intermediate r] of 
1.0 < r] < 2.8. 



4. Summary 

We present STAR's measurements of third harmonic flow (^3) of charged particles from Au+Au 
collisions at ^/snn = 200 GeV as a function of pseudorapidity, transverse momentum, and 
centrality. We report result rom a two-particle method for particle pairs with an gap or fit 
with a wide Gaussian in pseudorapidity separation, as well as from the standard event-plane 
method with an event plane near midrapidity or at forward rapidity. Short-range correlations are 
eliminated either by an gap or by eliminating the narrow Gaussian in pseudorapidity separation. 
The measured values of V3 continuously decrease as the mean pseudorapidity separation((Ar/)) 
of the particles increase within the range observable by STAR. It is not known whether this 
decrease is due to a decrease in non-flow correlations or due to a decrease in fluctuations. We 
observe that ^3 increases with transverse momentum before it levels-off above 3 GeV/c, similar 
to the case of elliptic flow. Our results are mostly described by hydrodynamic models with small 
viscosity with MC Glauber initial conditions. 
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Figure 8. V2 (top) and fs (bottom) for Au+Au collisions at ^snn = 200 GeV in 0-5% and 
30-40% central collisions as a function of transverse momentum, compared with ideal [10] and 
viscous [12] hydro, AMPT transport [13], and Parton-Hadron String Dynamics [20j models. The 
STAR V2 values (top) come from Ref. [3]. 
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